Paramecium swimming in capillary tube 
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Swimming organisms in their natural habitat navigate through a wide array of geometries and 
chemical environments. Interaction with the boundaries is ubiquitous and can significantly modify 
the swimming characteristics of the organism as observed under ideal conditions. We study the 
dynamics of ciliary locomotion in Paramecium multimicronucleatum and observe the effect of the 
soHd boundaries on the velocities in the near field of the organism. Experimental observations show 
that Paramecium executes helical trajectories that slowly transition to straight line motion as the 
diameter of the capillary tubes decrease. Theoretically this system is modeled as an undulating 
cylinder with pressure gradient and compared with experiments; showing that such considerations 
are necessary for modeling finite sized organisms in the restrictive geometries. 



I. INTRODUCTION 

Microorganisms use a variety of propulsion mecha- 
nisms to swim around in their habitat for predator eva- 
sion [I or locomotion towards favorable gradients [5]. 
The locomotory behavior and performance is usually con- 
trolled by chemical or hydrodynamic cues arising due to 
interaction with the local environment. The behavior of 
organisms can be studied considering them as a single en- 
tity or in a group. Collective motion of organisms have 
revealed vivid characteristics; for example sperm swim- 
ming near surfaces execute circular trajectories and ag- 
gregation of organisms near the surfaces is an interesting 
phenomena f3^. In addition, swimming velocity and di- 
rection, tumbling probability and turn angle in capillaries 
have been characterized [H [S] . Motility of cells and the 
morphological changes due to restrictive geometries are 
active area of interest [HI [7] . In all of these systems organ- 
isms often interact with the nearby surfaces of organisms 
or the boundaries (geometrical constraints) imposed by 
the nature on the motion; which causes the organism to 
exhibit varied swimming characteristics as compared to 
its motion in ideal infinite fluid medium. 

In many eukaryotic microorganisms coordinated mo- 
tion of cilia helps in propagating metachronal waves; 
which propel the organism[5]. Millions of ciliary hairs 
in mammals help in mucus transport and also function 
as sensory organelles that help in maintaining balance^. 
The ciliary beat is an interesting physical phenomena 
and has been studied extensively from biological point of 
view 1 101 . For example, experiments have been conducted 
to gain insights into the electro-physiological changes in 
the ciliary beat [TT] and taxis of ciliates under various 
conditions have been categorized [T2] . The effect of high 
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viscosities [13] on the locomotory traits of Paramecium 
especially with regards to changes in wave velocities, am- 
plitude of beat of the cilia, wavelength of the ciliary 
beat have been extensively documented. Experiments 
on swimming of ciliates in vertically aligned tapered glass 
tubes [H] have revealed interesting locomotory traits. 

A 2D wavy sheet that can be used to model simpli- 
fied swimming motion in many micro-organisms provided 
the theoretical framework for modeling microorganisms, 
by considering small amplitude expansions of the prop- 
agating wave [15]. Further studies into the swimming 
of 2D sheet surrounded by the planar boundaries re- 
vealed propulsive advantages for specific beat patterns 
[TB] . Phase locking in wavy sheets have been studied 
theoretically [17] and fiow patterns in the near-field of 
cilia have been investigated numerically [T5]. Thousands 
of cilia in Paramecium beat just out of phase to propa- 
gate waves in fluid and hence can either be modeled as 
infinitely long cylinder T9] or sphere ^20] with surface un- 
dulations or as a spheroid with slip velocity ^21j. These 
models were also used to validate swimming patterns in 
variety of ciliates [52] and to develop a boundary layer 
theory for predicting the near and far field velocities of 
ciliary micro-organisms in the unbounded fiuid |23] . 

Previous experiments have involved measuring the av- 
erage velocity; however various trajectories executed by 
Paramecium have not been considered. Theoretical stud- 
ies involving swimming Paramecium have mostly fo- 
cussed on infinite models without consideration of the 
boundary effects. We present a combined experimental 
and theoretical approach to reveal the locomotive pat- 
terns of Paramecium multimicronucleatum and rational- 
ize the pressure gradient effects, on the swimming behav- 
ior in confined spaces. In section|ll| we explain the exper- 
imental methods of introducing the organisms in confined 
geometries like the capillary tubes and the techniques 
used to visualize their motion. A theoretical model with 
a pressure gradient is developed to understand the effect 
of boundaries in section |III[ In section |IV[ we compare 
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the predicted swimming velocity with the experiments 
and discuss other important parameters that might affect 
locomotion of microorganisms close to the boundaries in 
section |Vl 



II. EXPERIMENTS 
A. EXPERIMENTAL METHODS 

Paramecium multimicronucleatum is a single-celled eu- 
karyote commonly found in warmer regions of the fresh- 
water ponds [8 . Cultures were reared in a double wheat 
medium and subcultures were placed every 11 days when 
they reach their peak population. Paramecium at the be- 
ginning of the their exponential growth curve were used 
for the experiments. The cultures were centrifuged and 
washed twice in a Buffer solution consisting 9 mM CaCl2, 
3 mM KCl, 5 mM Tris-HCl (pH 7.2) to remove the debris 
and were allowed to equilibrate for 30 min. The equili- 
brated cells were then observed under the Leica DMI 
3000 microscope at 5X, 20X and 40X with bright field or 
DIG optics and their motion was recorded using Redlake 
MotionXtra N3 camera. 

The ciliary coordination in ciliates is often controlled 
by a complex collection of external cues that causes the 
organism to change the frequency or other parameters 
of wave propagation. However the difference in direc- 
tions of the propagating metachronal wave and swim- 
ming direction causes the organism to move in helical 
path. Broadly the locomotory gait can be classified 
as forward(anterior) or backward(posterior) swimming; 
with the forward swimming exhibiting different helical 
modes. The anterior swimming direction and wave prop- 
agation direction are separated by 135° in right handed 
helix swimmers and by 225° in left handed helical swim- 
mers [12] . The ciliary reversal modes of locomotion are 
characterized by little or no helical motion and a coasting 
motion with lower swimming speed. The helical modes 
of swimming in contrast to the ciliary reversal modes 
allows us to better characterize the change in the loco- 
motive pattern and hence will be used to study the effect 
of boundaries. 

To investigate the fiow-field around the organism sus- 
pensions of Polystyrene microspheres (5 /zm diameter, 
Thermo Scientific) prepared in EDC solution were intro- 
duced into the cultures. Small volume of the cultures 5 /zl 
were then placed on the glass slide which created a very 
thin film and allowed us to visualize the 2D flow field 
around the Paramecium. Figure 1(a) shows a Parame- 
cium swimming in a thin film of liquid. We can see there 
are two strong vortices forming on the lateral sides of the 
organism; showing the strong tangential velocity of cilia 
on the far field of the organism. 

The effect of confined geometries on the ciliary dynam- 
ics is examined by introducing the organisms in capillary 
tubes. Tubes required for this purpose are manufactured 
by attaching a dead weight to the end of the borosilicate 




FIG. 1: (a) A Paramecium swimming in a thin layer of fluid, 
(b) Streaklines of particles in flows generated by the cilia 

glass pipettes and by heating their tips. By controlling 
the value of the dead weight and the intensity of the 
applied heat different diameters of capillary tubes rang- 
ing from 90^250 fj,m were manufactured. Some commer- 
cially available tubes with specific diameters(d=100,150, 
200 /im) made of borosilicate glass were ordered from Vit- 
rotubes. The equilibrated cultures were then transferred 
to the extruded glass pipettes where they got pulled into 
the small constant cross section of the tube due to capil- 
lary forces. 

B. EXPERIMENTAL OBSERVATIONS 

Paramecium were found to have a long and slender 
structure with typical lengths around 212±14 and the 
width about 57±5 fim (shown in Fig. [l|a)). The velocity 
of these micro-swimmers in unbounded fluid was found 
to be 1064±83 fim/s. We then performed experiments 
with capillary tubes of different diameters. Paramecium 
swimming in buffer (isotonic solution) were put into the 
capillary tubes which caused them to be confined in small 
circular geometry. A generic code written in MATLAB 
was used to track the motion of these organisms and the 
velocities were computed. 

In order to measure the vital parameters for swimming; 
we captured the cilia motion with the high speed camera 
at 300 fps, which allowed us to visualize the metachronal 
wave propagation over the organism. Each cilium was 
found to be 10~12 /im in length and 0.2 fim in diameter 
and beats slightly out of phase compared to the nearby 
cilium, thereby causing a traveling wave to pass over the 
surface of the organism. The typical wavelengths of the 
metachronal waves measured from our experiments were 
27 /im, half peak to peak amplitude 4.2 /zm and the fre- 
quency of the beat being around 30 Hz We assume 
that the cilia are so closely packed that the fluid does 
not penetrate the material wave, thereby allowing no slip 
boundary condition to hold good. 

Imaging of round capillaries under microscope caused 
optical distortions which leads to recording of altered am- 
plitudes and velocities. We directly took images of the 
cross-section of various capillaries to get a relation be- 
tween the true and the observed inner diameters of the 
tubes, which was further used to correct the observed 



3 



amplitude and velocities of the organisms. 




FIG. 2: The metachronal waves propagated by the ciha of the 
Paramecium over a period of 33.3 ms. The sweep of the cilia 
and the direction of propagation of the wave can be visuahzed. 
Arrows follow the peak amplitude of the waves. 

In tubes of extremely small diameter the swimming ve- 
locity of the organism was very low with almost a straight 
line motion with variable rotation rates. Whereas in 
tubes of larger diameter the organism was seen to move 
in a helical path instead of straight line motion. We ap- 
plied a correction factor for the path of the swimming 
organisms, as the image analysis only revealed the 2D 
projection of the helix. 

It is observed that the Paramecium swims slowly as the 
tube diameter is decreased. This can be attributed to the 
increased drag felt by the organism due to the proximity 
of the boundaries. For the capillary tubes whose diame- 
ter were very close to the diameter of the Paramecium^ 
the swimming velocity was close to zero. 

In tubes of smaller diameter {R/c ^ 1.5) we observed 
that a backward (posterior) swimming Paramecium ex- 
ecuted a helical swimming trajectory with small ampli- 
tude wavelengths. Such swimming gait have not been 
reported before [HI . In this range of tube diameters 
and for the forward swimming Paramecium, little or no 
helical motion of the organism is observed. 



III. THEORETICAL MODEL 

The governing equations for very low Reynolds num- 
ber flows {Re ^ 1) are the Stokes equations: Vp — 
/iV^u, V • u = 0. Since these organisms have a large 
length to diameter ratio; they can be effectively modeled 
as cylinders. As seen in the Figure 2 the cilia create syn- 
chronized motion to produce metachronal waves; to an 
observer this seems like a material wave propagating on 
the surface of the organism. Under the assumption that 



no fluid penetrates the wave of cilia tips since they are 
densely packed, the concept of envelop model can be ap- 
plied for the physical system. The boundary in this case 
is the circular capillary tube; thus the problem reduces 
to modeling a cylinder with a wavy surface swimming 
inside a cylindrical geometry. 

Figure |4] shows the schematic of the organism swim- 
ming the tube with a velocity Uswim- Paramecium has 
previously been modeled as an infinitely long cylinder 
so that there is no pressure gradient at the front and 
the back while swimming in the unbounded fluid. Sim- 
ilar infinite models can be developed for a Paramecium 
swimming inside a tube [19]. However, due to the pres- 
ence of the confined spaces the finite length Paramecium 
experiences a pressure gradient at its ends that influences 
its locomotion within tubes. Considering the propagating 
wave to be of frequency w, wavelength A, radial amplitude 
6, transverse amplitude a, and wave number fc = ^; any 
material point S on the undulating surface can be written 
as: 

rg = c + bsai{kz — ujt), Zs — z + acos{kz — cot + (jj) (1) 
where (j) is the phase difference. 

We work in the frame in which the organism is swim- 
ming with a velocity Uswim- In the small amplitude limit, 
the boundary conditions of the surface of the organism 
at r = c and on the solid walls at r = i? can be written 
as: 

Ur\r=c — — fewCOs(fcz — Ujt),Uz\r=c — flW sin(fcz — Ujt + (f)) 
Ur\r=R — 0,Uz\r=R = UswiTU (2) 

Thus the swimming problem can thus be envisioned as 
sum of pressure driven flow and shear flow in the nar- 
row annulus that surrounds the Paramecium. We seek 
a solution in terms of the streamfunction -0 such that 
■0 = ip(p) -\- where i/j^pHs the streamfunction corre- 
sponding to pressure driven flow and '0*^^Mue to the the 
shear flow. Using cylindrical co-ordinates and axisym- 
metric potential theory we can write the velocity com- 
ponents to be Ur = —-f^jUz = Solving for the 
pressure driven flow with stationary boundaries we can 
get the streamfunction to be: 



(p) 



1 dp 
4fi dz 



R^- 
In^ 



2 2 



(3) 



For the shear flow part we can substitute the ve- 
locity components in the Stokes equation and by tak- 
ing curl we end up with a equation of the form 



r dr 



Rewriting coordinate systems as tj = kr and C = fcz — 
Lut and using separation of variables in r/ and C. we can 
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(a) R/c=2.63 (b) R/c= 1.67 (c) R/c=0.90 



0.5 mm 



0.4 sec 
► 



0.8 sec H^^^^MHMMBIKl^^^^^^^V 

^ • • • - 

► - • • 



FIG. 3: Swimming of Paramectum in tubes of different diameters. Fig. (a) shows swimming in large tube (R/c =2.63) where 
the trajectory of the motion is heUcal. In swimming inside tubes of intermediate diameters (R/c=1.67 ) small wavelength 
helices are seen. In very small tubes(R/c = 0.9) as shown in Fig. (c) Paramecium swims in a straight line. 



Undulating wave 




FIG. 4: Schematic of a wavy Fig. (a) infinite cylinder or 
Fig. (b) finite cylinder swimming inside the cylindrical 
tube(boundaries) . It is assumed that the organism is swim- 
ming in the negative z direction with a velocity Uswim- By 
transforming the coordinate system to one, where the or- 
ganism is swimming with velocity Uswim we get the current 
schematic with the descriptions of the velocities on the walls 
of capillary. 



with Equation [2j For calculations of the zeroth and first 
order the Un, Pn,Jn,Sn does not have any contribution 
and hence we are left to determine four unknowns for 
the problem. 

Coupled with the pressure driven flow this gives 
us a set of five equations and six unknowns 
Ai, Bi,Ci, Di,Uswimj Also from the flux condition 
of the swimming Paramecium we have Q — Uswim'^^ — 
J J UzrdrdO. This gives us a condition between Uswim 
and . Thus we have a set of five equations and five un- 
knowns which can be solved to find out the constants and 
the swimming velocity Uswim of the organism. The ex- 
pression for swimming velocity of the Paramecium comes 
out to be: 



a^kuj abkric cos dp 

2 ^ s^i ac 



abk cos I 



(5) 



obtain a streamfunction solution of the form; 



(^) - 



^F„sin(C)+ VG„cos(C) (4) 



i=0 



=0 



where, 

Gn = V [anKi{nri) 



The above expression shows that the swimming veloc- 
ity is directly dependent on (afc)^ and also the pressure 
gradient terms. For the infinite boundaries case with no 
pressure gradient the above equation reduces to the fol- 
lowing expression: 



and An,Bn, 



BnriK2{nri) + Cnhini]) + Dnrjhinfii)] 
l3nr]K2{nr]) + 7„/i(nr;) -|- SnVhinrj)] 
Pmln,6n are the constants to 
be determined from the boundary conditions and / and 
K are modified Bessel functions of the first and second 
kind. 

We seek perturbation expansions of the velocities Uj. — 



.fc2^ and 



^'^'^drj derived from the stream func- 



tion to the 0^{bk) and 0^{bk) and compare the velocities 



Ublak 



Lok 
~2 



{K',^Kl)b^ 



Kl 



K0K2 

which is same as in 1191. 



2Klab COS0 



(6) 



IV. RESULTS 

The motivation of the study was to rationalize the be- 
havior of the organisms in close proximity to the bound- 
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aries. As the Paramecium swims inside the tube it traces 
out a hehcal path with the length of the body being 
ahgned in the swimming direction. Figure [5] shows the 
variation of the amphtude of the organism as it swims 
in capillary tube of different diameters. Theoretically we 
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FIG. 5: Plot of the non-dimensional amplitude of the helical 
path vs the non-dimensional radius of the tube. The plot 
shows a linear variation which increases with the radius of 
the tube. 

considered the small amplitude expansions for the waves 
and hence the expression derived for the swimming ve- 
locity would only be valid for smaller capillary tubes in 
which there is very less off axis movement. In larger tubes 
the organism does remain very close to the glass surfaces; 
but there is asymmetry in the boundary conditions that 
need to be applied. From the plot above and the con- 
straints on the amplitude of swimming of Paramecium 
we can conclude that for A/ {2c) < 1 our experimental 
and theoretical results would remain valid. This gives us 
corresponding value of the non-dimensional radius of the 
capillary tube R/c <2. 

Figure [6] shows the variation of velocity with the ra- 
dius of the capillary tube. The dotted line shows the 
swimming velocity when no pressure gradient effects are 
considered and the solution converges very quickly to the 
case for which the boundaries are at infinity. In the ex- 
periments we observe a slowly increasing trend of veloc- 
ity. For a finite size Paramecium swimming inside a re- 
strictive geometry there is a finite pressure gradient at 
the front and the back. The plot of the swimming veloc- 
ity considering finite pressure gradient is shown by the 
solid line. It shows a slowly increasing trend and finally 
converges to the infinite boundary and no pressure gra- 
dient case for very large R/c values. 

It can be seen that for R/c < 2 both the experi- 
ments and theory show an increasing trend and predic- 
tions match quite well. For the larger R/c we see that 
the velocities are much larger as opposed to that pre- 
dicted value for the infinite boundary case. From the 
experiments we see that the organism follows the surface 
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FIG. 6: Plot of normalized swimming velocity vs. the no- 
dimensional radius of the tube. Two different cases are com- 
pared, one with the dotted line shows the model of an infinite 
ciliate swimming in capillary tube. The solid line shows the 
velocity variation of the finite sized Paramecium inside con- 
fined geometries. Dots represent the experimental result and 
compares well with the finite pressure gradient model due to 
the restrictive environment. 



of the glass capillary as it swims in helical path in tubes 
of larger diameter. This situation can be thought of as 
Paramecium swimming close to a single solid wall and 
hence explains the observation of larger swimming ve- 
locities for tubes of larger diameter; where it constantly 
swims close to the wall. This reflects the large errors 
corresponding to the swimming case in larger tubes. 
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FIG. 7: Plot of the non-dimensional wavelength vs. non- 
dimensional radius of the tube. No apparent trend is seen 
with increasing radius of the tube. 

Figure [7] shows the variation of wavelengths of the he- 
lical path while they are swimming in the tube. The 
wavelengths also show a linearly increasing trend for the 
limit R/c < 2. In this limit the helices of the swimming 
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are well defined and periodic wavelengths are observed. 
While swimming in tubes of larger diameter(i?/c > 2) 
it was observed many times that the organism did not 
execute a full helix. 



V. DISCUSSION 

Wc investigated the locomotion of ciliary organisms in 
confined geometries. As the boundaries close on the or- 
ganism, more viscous effects are felt by the cilia. The 
swimming velocity of Paramecium decreases due to the 
effect of close boundaries. Also in such confined spaces a 
finite sized organism feels a pressure gradient across the 
ends. This pressure gradient affects the swimming ve- 
locity and needs to be considered while modeling similar 
self propelling objects in restrictive geometries. Our ex- 
periments confirm this observation which show a slowly 
increasing trend of velocities. 

Many interesting questions arise from this study espe- 
cially about the amplitude of the waves propagated by 
the Paramecium. In a restrictive channel the beat of the 
waves is limited by the dimensions of channel and the size 
of organism. Due to the proximity, do the amplitudes of 
the beat change to provide a locomotory advantage? 



The helical path traced out by the Paramecium with 
the anterior portion of the body aligned towards the local 
swimming direction is also an interesting locmotory trait. 
It was also observed that as the radii of the capillary tube 
increased the radius of the swimming helix also increased. 
It might be possible that being close to the curved surface 
provides a propulsive advantage. It was also found that 
anterior swimming Paramecium can execute well defined 
helical paths when put inside capillary tubes of certain 
diameter. 

The study revealed the interesting locomotory traits 
in presence of solid wall. Future work would involve in- 
vestigating the hydrodynamic effects of different textured 
boundaries on the swimming characteristics of the organ- 
ism. 
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